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ABSTRACT 

A sample of faint, V magnitude selected, galaxy pairs, having physical separations less than 
approximately 20h~^ kpc, is used to examine the rise in the merger rate with redshift and the sta- 
tistical relations between close pairs and the field galaxy population. Redshifts have been obtained 
for 14 galaxies {V < 22.5) that are in close {9 < 6") pairs, along with a comparison sample of 38 
field galaxies. Two color photometry is available for about 1000 galaxies in the same fields. The 
average redshift of the V < 22.5 field population is 0.36, statistically equal to the average redshift 
of 0.42 for the pairs. The similarity of the two redshift distributions, Az < 0.1, limits any differen- 
tial luminosity enhancement of close pairs to less than half a magnitude. The pairs are somewhat 
bluer than the field and have nearly twice the average [O II] detection rate of the field, but the 
differences are not statistically significant. The field population has an angular correlation at sepa- 
rations of 9 <6" higher than the inward extrapolation of uj{9) oc 0~°'*, which may be a population 
of "companions" not present at the current epoch, or, luminosity enhancement of intrinsically faint 
galaxies in pairs. Physical pairs comprise about 7% of the faint galaxies in our survey fields. The 
same physical separation applied to local galaxies finds only 2.6% in pairs. If the rise in close low 
relatively velocity pairs with redshift is parameterized as (1 + z)™", then m = 2.9 it 0.8. If all pairs 
at low velocities and r < 20h~^ kpc merge, then the average galaxy mass would be 32% smaller at 
z = 0.4 than locally. 



1. INTRODUCTION 

Pairs of galaxies at separations of less than 10-20/i~^ kpc comprise a small fraction, about 
3%, of the present day galaxy population. Many close pairs of galaxies (e. g. Arp and Madore 1986) 
are disturbed, sometimes grossly, by gravitational interaction with a neighboring galaxy (Toomre 
1977, Larson and Tinsley 1978). The distorted morphologies can be understood if the infalling 
galaxies are on orbits close to parabolic, which, in the presence of dynamical friction, implies that 
most close pairs will merge to produce a single galaxy within 10^ years (e. g. Toomre 1977, Barnes 
1988). The impact of merging on galaxy numbers and masses is a small correction, of order 10%, if 
it occurred at a uniform rate in the past; however, a number of straightforward arguments predict 
a rapid rise of merging with redshift within any gravitational instability theory for the formation 
of galaxies. Toomre (1977) shows that if the distribution of marginal pairwise binding energies 
is a constant, then the rate of binary mergers increases into the past as t~^'^ , corresponding to 
(1 + z)^' ^ in an $7 = 1 Universe. Somewhat more general considerations, using the Press-Schechter 
(1974) formalism and including a decoupling from the merging hierarchy, indicate that the rate 
of merging is proportional to (1 + z)"*, where ra ~ 4.5rj'^'^^ at low redshift (Carlberg 1990) for a 
fixed mass spectrum. Extrapolation of present day merger rates using functions with such a strong 
redshift dependence can have a dramatic effect on the galaxy population; such considerations imply 
that L* galaxies might be relatively rare beyond z = 1 (Broadhurst, Ellis and Shanks 1988; Rocca- 
Volmerange and Guiderdoni 1990; Carlberg and Chariot 1990; Cowie, Songaila, and Hu 1991; 
Broadhurst, Ellis, and Glazebrook 1992; Carlberg 1992), although this conclusion is very sensitive 
to the normalization and redshift dependence of the merger rate. 

Zepf and Koo (1989) undertook a study of faint galaxy pairs based on KPNO 4m plates. 
They found 39 paired galaxies with separations less than 4.5", drawn from a total population of 
1055 galaxies with mj < 22. A comparison with low redshift pair samples led them to conclude 
that the merger rate rises as (1 + z)"^ , under the assumption that the pair redshift distribution 
is similar to that observed in the field; furthermore, they noted that many of the pairs showed 
evidence of distortions that might arise from tidal effects. A recent study of close pairs using HST 
data (Burkey etal. 1994) finds that the merger rate depends on redshift as (1 + z)^'^, in agreement 
with the above. Indirect evidence that the rate of merging does indeed rise rapidly with redshift 
can also be inferred from a variety of galaxy phenomena believed to be related to mergers - for 
example, quasars (Boyle etal. 1988), IRAS galaxies, (Lonsdale etal. 1990), and the Butcher-Oemler 
effect (Dressier and Gunn 1983, Butcher and Oemler 1984, Lavery and Henry 1988, Lavery, Henry, 
and McClure 1992). 

The redshift distribution of a magnitude-limited survey of pairs of galaxies could be quite 
different from that of the field, for at least three reasons. If paired galaxies are, on average, 
more luminous than field galaxies (perhaps as a result of induced star formation) then their mean 
redshifts will be higher. However, a magnitude-limited survey of pairs of galaxies that are about to 
merge may have somewhat lower mean masses than field galaxies (if pair luminosity enhancement 
is primarily concentrated in intrinsically low luminosity galaxies). On the other hand, if pairs have 
the same luminosities as field galaxies, then selecting pairs at a fixed angular separation favors 
pairs at low redshift, because S,i9dA{z)), the real space correlation function, decreases along the 
line of sight. (Here dAiz) is the angular diameter distance.) For the limiting case of no luminosity 



evolution, the numbers of physical pairs per unit volume, np{z), are distributed (for an $7 = 1 
Universe) as 

Upiz) oc (1 + z)-^-'+^[l - (1 + zy-^j'-'n^z), (1) 

where n{z) is the redshift distribution of single galaxies, the correlation function is modeled to evolve 
as {1 + z)~^~^ (e. g. Efstathiou etal. 1991), and the real space correlation function ^ oc r~'^. For the 
relatively shallow depths of interest here, z ~ 0.3, the average redshift of the pair distribution is 
expected to be about 10% smaller than that of the field. We conclude that the average redshift of 
pairs, compared to that of the field, may be an interesting constraint on any differential luminosity 
evolution of the pairs. (Nevertheless, the effect of equation (1) could potentially be a 50% reduction 
in mean redshift for samples that have a significant fraction of their population spread beyond 
z>l.) 

The main goals of this paper are to compare the redshift distribution of close pairs and 
the field, and to estimate the rise in the merger rate with redshift using a CCD selected sample 
of pairs both with and without redshifts. The next section discusses our selection procedure and 
observations. Section 3 contains our observational results (redshifts, colors, and angular correlation 
function) . The implications of these data for the evolution of galaxies in mass and luminosity are 
discussed in S4. 



2. PAIR SELECTION AND OBSERVATIONS 

There are several straightforward considerations in constructing a sample of galaxy pairs for 
a small redshift survey. The sky density of physical pairs with separations less than some small 
angle 6 is Np = -/Vo(/)/o (^ + i^{0))2'K9d6 at flux level /. The angular correlation function, uj{9), 
can be accurately represented as a power law with slope near 1 — 7 ~ —0.8. The amplitude scales 
with depth in a approximately Euclidean manner (Maddox et al. 1990, Pritchet and Infante 1992) 
- i. e. as uj{d) oc 0~t+i/')'/2. xhe sky density of faint galaxies brighter than magnitude m rises as 
logn(< m) oc 0.45m, or no(> /) oc /~^'^^^. In the regime where uj{0) > 1, the surface density of 
physical pairs at fixed angular separation 9 is proportional to nQ(> f)uj{6), or Up oc f~^'^^ ~ i. e. the 
number of pairs brighter than m obeys log A'p oc 0.54?n,. Thus the close pair count rises somewhat 
faster with magnitude than the number of galaxies - a fact which has an interesting implication 
for the optimal magnitude at which to observe faint physical pairs. For some fixed total observing 
time T with a multiobject spectrograph, the optimal exposure time t on each field is clearly that 
which maximizes the rate R = Up/t at which useful pair spectra are acquired. Assuming that all 
spectra above some limiting continuum signal-to-noise ratio have measurable redshifts, then one 
finds that, in the read noise limit, the rate R is given by i? oc t*^"*^^, whereas in the sky noise limit, 
R oc t"*^-^^. Therefore the total number of pairs will be optimized for pairs selected brighter than 
the magnitude at which the spectra begin to be sky-noise dominated {V ~ 22 in our case). 

The field galaxies at y ~ 22 have a median redshift of z ~ 0.35 (CoUess etal. 1990, Lilly etal. 
1991), sufficient for some redshift leverage on the merger rate-redshift relation. At this redshift, 
the field galaxy population is significantly different than that observed locally: it is bluer, and the 
luminosity density is about 3 times higher (Colless etal. 1990, Lilly etal. 1991, Eales 1993) (even 
though the characteristic luminosity L* is about the same as observed locally). 



For our purposes pairs should ideally be defined as those objects that are doomed to merge; 
however, for projected data in a continuous clustering hierarchy the definition of a pair is somewhat 
arbitrary. A reasonable operational definition is to select pair candidates at angles such that 
uj{9) ^ 1 {i. e. > 50% probability that the pair is physical and not projected). At y = 22.5 the 
amplitude of the angular correlation function is loguj{l°) ~ —2.5, or uj{l") = 2 (Pritchet and 
Infante 1992), implying that uj{9) > 1 for 9 < 3" (for uj{9) oc 9~^'^). The pair density on the 
sky is quite low: the expected number of paired galaxies with 9 < 6" (calculated by extrapolating 
the Pritchet and Infante 1992 correlation function) is about 1.5 pairs per 6x5 arcminute field 
available for spectroscopy. To increase our chances of obtaining pair spectra, pairs were preselected 
from the same plates as used for the "CFHT North Galactic Pole Survey" (Infante and Pritchet 
1992). Because of concerns regarding the separation of close pairs in the original APM image mode 
data, we ran the Kron (1980) galaxy finding and photometry algorithms on scans of the central 
regions of the plates. To obtain a reasonable number of pairs we selected all paired galaxies with 
separations 9 < 8" (rather than the 3" criterion discussed above). Based on the Pritchet and 
Infante (1992) correlation function, one might anticipate that ~ 4 of all pairs with separations of 
8" and V magnitudes < 22.5 will be physical. For V < 21.5 this fraction will rise to ^ 60%, and to 
~ 80% for pairs chosen with a separation in the range 2-8". 

Chosen pairs were all brighter than J = 22.5 in photometry done on the APM scan data 
from Infante and Pritchet (1992). We then selected fields to (i) maximize the number of pairs (at 
least 2 pairs per field, never more than 3), and (ii) preselect (where possible) pairs with similar 
orientation for spectroscopy. Note that this selection procedure will bias upwards the number of 
small separation pairs. Our original selection procedure (magnitude, separation, and visual check) 
identified 110 pairs in an area equivalent to 80 spectroscopy fields (i. e. 1.4 pairs per field). On the 
basis of the average numbers of pairs selected (--^2.5 per field) and expected (~1.4 per field), our 
procedure should result in a boost of approximately 60% in the number of close pairs over that 
expected in purely random fields. However, we show in §3 that the presence of these preselected 
pairs in our fields actually has only a small effect on our computed correlation functions (because 
the correlations of close pairs turn out to be much stronger than expected). 

Using the above considerations to choose fields, observations were acquired at the Canada- 
France-Hawaii Telescope in 1992 March using the Marlin multi-object imaging spectrograph in 
conjunction with a thick, front-side illuminated, Ford Aerospace/SAIC 1024^ CCD (read noise 6.4 
e~ pixel"^; peak quantum efficiency qmax — 40% at 6500A; q = \qmax at 5000A and 8000A). 
Imaging observations were acquired of seven of our preselected fields, each approximately 6' x 6' in 
area, through Johnson V and Kron-Cousins / filters (exposure times of 15 min in each filter). 

The imaging data were used to design a slitlet mask for each field, using software available 
at the telescope. First, preselected pairs (see above) were identified on the images, and as many 
pair galaxies as possible were assigned a slit width (1.5") and length (nominally 15"). (It was not 
possible to obtain spectra of all galaxies in pairs, since >50% of them had orientations that would 
have resulted in overlapping spectra.) Once the pairs had been assigned slits, the rest of each field 
was filled with field objects, many of which were below our expected spectroscopic completeness 
limit oiV ~ 22. The final data file containing slit coordinates was sent to a computer-controlled 
laser punching machine ("LAMA"), which produced the final slitlet mask. Spectroscopy was then 
obtained over a useful wavelength range AA5000-8000, with a dispersion of 7.7A pix~^ and slits 



1.5" wide and usually 10"-15" long. The total spectroscopic integration time per field was 2 hours, 
spread over 4 individual 30 min exposures. Short (60 s) direct images were taken through the mask 
(at the beginning and in the middle of each 2 hour sequence of exposures), to check the alignment 
of objects in the masks, and hence test for flexure. Arcs were acquired at the beginning and end 
of each 2 hour integration on a field. 

Spectra were extracted using the IRAF^ apextract package. Redshifts were obtained by 
identifying two or more of [Oil] A3727, 4000A break, H/3 A4861, [OIII] AA4957,5007, Mg b A5180, 
and in some cases Ha A6563. The continuum shape also played a role in checking redshifts obtained 
from features and breaks. A quality index Q, ranging from 1 (highest quality) to 6, was assigned to 
each redshift. Most of the analysis below refers to redshifts with Q < 3, although including (5 = 4 
redshifts changes none of our conclusions. 

The imaging fields were later analyzed using photometry and image classification loosely 
based on the Kron (1980) algorithm. At y < 22.5 (and / < 23) our survey contains 1062 objects, 
of which about 60% are galaxies, in an area of about 9 x 10^ square arcseconds (some area was 
lost due to internal reflections from bright stars). Internal photometric errors are about 0.1 mag 
at V=22.5 and 1=22. The image quality is good, 0.8", in the field center, but degrades to 1.5" in 
the field corners. 



3. OBSERVATIONAL RESULTS 

The average properties of the pairs (chosen to have 9 < 6") and field sample are given in 
Table 1. The second column gives the number of objects placed in slits (excluding stars and 
bad extractions), and the third column gives the number of objects, Nz, for which we were able to 
obtain redshifts. The pair redshift sample is 8/14=57% complete, and the field is about 28/63=44% 
complete (for V < 21.5). There are at least two reasons for this low completeness, (i) The CCD 
that we used had poor sensitivity in the blue, with essentially no response below --^ 4500A, and only 
~ 4 of peak response at 4800A. This made the 4000A break and [O II] A3727 impossible to detect 
for z < 0.1 objects, and very difficult at z < 0.2. Interpolating between the results of Colless etal. 
(1990) and Broadhurst etal. (1988), we find that approximately ^ of the galaxies in our V < 21.5 
sample should lie at redshifts z < 0.2, compared with 5/38=13% actually observed. Our overall 
completeness for objects brighter than V = 21.5 and with z > 0.2 therefore rises to ~ 60%. (ii) 
The information in Table 1 is a fairly conservative cut of objects with secure redshifts. Including 
the next quality level of redshifts increases the completeness to 61% at F < 21.5, or ^^ 80% for 
objects with z > 0.2. (We note in passing that such a change in redshift quality has hardly any 
effect on the mean redshift, mean color, or mean line fluxes quoted in Table 1.) 

The fourth column of Table 1 gives the average redshifts for V < 21.5 and V < 22.5; these are 
a little higher than the values of 0.25 to 0.35 expected for the two magnitude ranges (Colless etal. 
1990). Apart from one extragalactic HII region, the lowest redshift detected is 0.137, a consequence 



^ IRAF is distributed by the National Optical Astronomical Observatories, which is operated 
by the Association of Universities for Research in Astronomy, Inc., under contract to the National 
Science Foundation. 



Table 1: Average Properties of Paired and Field Galaxies 



Sample 


NsHt 


N. 


(z) 


(V) 


{V-I) 


(L(0 II)) 
[h ^ erg s ^] 


V < 21.5 pairs 

V < 21.5 field 

y < 22.5 pairs 

V < 22.5 field 


14 
63 

24 
97 


8 
28 

14 
38 


0.28±0.03 
0.36±0.03 

0.42±0.05 
0.37±0.03 


20.5 
20.7 

21.1 
21.1 


1.19±0.13 
1.34±0.09 

1.41±0.15 
1.33±0.09 


2.4 ± 1.3 X lO-^o 
1.8 ±1.1 X 10"^° 

4.6 ± 1.9 X 10"^° 
2.3 ± 1.4 X 10^0 



Note: all errors are errors in the mean. 



of the poor sensitivity of our CCD below 4500A, as discussed above. The dispersion about the mean 
redshift is about 0.16 for both field and pair samples. 

3.1 Redshift Distribution 

The most significant result arising from the redshift data is that the average redshifts of the 
pair and field samples are the same to within Az/z < 0.1. This requires that both populations 
have similar characteristic luminosities: the intrinsic luminosity of the two samples is identical to 
within 40%. (Formally the V < 21.5 sample places an even more stringent constraint on luminosity 
evolution of the pairs sample, because {z)pair < {z) field- Nevertheless, this result must be viewed 
with some caution because of the small number of redshifts for V < 21.5 pairs, and because of 
subclustering in these redshifts.) Overall, our redshifts place a powerful constraint on differential 
luminosity evolution between the two populations. (As noted in §1, the bias introduced into pair 
selection because of the correlation function is much smaller than this.) 

3.2 Colors and Luminosities 

Column 5 of Table 1 gives the average magnitude of the galaxies, and column 6 gives the 
mean color and the standard deviation of the mean. In the V < 21.5 sample the pairs are bluer 
than the field, but the difference is at the la level. (The colors have not been /c-corrected, but 
this should have little effect on the above result because the mean redshifts of the samples are the 
same.) Color differences of \A{V — /)| > 0.2 magnitudes are nevertheless ruled out. 

The average luminosity in the [O II] 3727A line, column 7, is calculated (assuming that non- 
detections have zero line luminosity) in an Hq = lOO/i, $7 = 1 cosmology. There is no significant 
difference in this quantity for the field and pair samples; this result is perhaps not too surprising, 
given the fact that fewer than 40% of the pairs and 20% of the field galaxies have detections. 
Perhaps the only difference that is notable in the [O II] A3727 properties of the survey is that 
detections are twice as frequent in the pair sample as in the field sample (5/14=36% for pairs 
versus 13/70=19% for field objects at V < 21.5; 9/24=38% for pairs versus 16/104=15% for field 
objects atV< 22.5). 

There are 7 pairs for which both members of the pair possess redshifts; 5 of these pairs 
belong to our "pre-selected" sample, whereas 2 were "serendipitous". (One "pair" is actually a 
triple.) (Since our sample of pairs was selected using the [very faint] cutoff magnitudes V < 22.5 
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and I < 23, and because of orientation effects, it is not surprising that many pairs possess only a 
single redshift.) Three of the pairs with full redshift information have Az < 0.002, one pair has 
Az = 0.031, and the other two are separated by Az > 0.2: that is, three pairs are physical, whereas 
four are projected. The redshifts of the physical pairs are 0.378, 0.412, and 0.586. Considering that 
pairs were selected on the basis that uj{6) > 0.5, it is entirely to be expected that of order half of 
the pairs should be physically associated. 

Table 2: Angular Correlations 





Vf < 21.5 


Vf < 22.5 


// < 21.5 


uj{< 6") 


4.80 ±1.19 


2.06 ±0.41 


1.07 ±0.29 


oj{&' - 24") 


0.25 ±0.14 


0.09 ± 0.06 


0.18 ±0.05 


uj{< 6")/lj(6 - 24") 


19.2 


24.0 


5.9 


N,{< 6") 


24 


56 


52 


Np{6" - 24") 


80 


306 


514 


iVp(<6")/iVp(6-24") 


0.30±0.07 


0.18±0.03 


0.10±0.02 



3.3 Angular Correlations 

Table 2 gives the data for the small separation end of the angular correlation function, derived 
from our imaging observations. The samples marked with a subscript / have been selected to be 
"fuzzy": that is, their Kron (1980) r_2 parameters have been preselected (using r_2 > 1.8 for this 
data) to remove (unclustered) stars. (Note that there may be some compact galaxies which are 
excluded by this criterion.) 

A detailed calculation shows that the 6 < 6" bin should have about 15 x fewer pairs than the 
larger 6 < ^ < 24" bin (assuming a power-law angular correlation function with power-law slope 
6 = 0.8 and uj{3") = 1, and that galaxies closer than 3" cannot be distinguished). (Note that we 
find no pairs with separations < 1.5".) The observed ratios of pairs are in fact ~3-5 times larger 
than this for the V samples (with the larger excess arising from the most reliable Vf < 21.5 sample). 
The correlation function, ij{< 6"), is expected to be 2.6 times larger than cj{6" — 24") (under the 
same assumptions as above). In fact the observed correlations in the V band are a factor of 7-9 
times larger than expected in the innermost bin (a factor of 7 for the V < 21.5 sample). Excess 
correlations for the // < 21.5 sample are smaller, but still statistically significant (approximately 
a 2(T excess in the pair ratio Np{< 6")/Np{6 — 24")). The lower excess correlations in the / band 
could be due to the higher mean redshift of this sample, or to the fact that the pairs are slightly 
bluer than the field (at least for the V < 21.5 sample). 

These excesses of close pairs (over what would be expected from an extrapolation of uj{9) at 
larger 9) are even more significant if 5 < 0.8 (e. g. Maddox etal. 1990), or if (as might be expected) 
there is incompleteness in the 9 < 6" bin. (At the mean redshift of the sample, 6" corresponds to a 



physical separation of 19h~^ kpc for fi = 1.) It is interesting to note that there is a weak indication 
of excess small scale correlation at similar separations at low redshifts (Davis and Peebles 1983), 
although the statistical significance of this result is unclear. 

The bias that would result from our field selection (§2) appears to fall far short of explaining 
the observed excess of close pairs. Only two of our preselected pairs survive the separation {A6 < 
6") and morphology (r_2 > 1.8) cuts that were used in Table 2; removing these two pairs would 
have no effect on the significance of the excess of correlation power at small separations. 

The images of all the close pairs were examined for signs of distortions and interactions. 
Although some of the galaxies have features that could be created in tidal distortions, few of the 
galaxies are strongly distorted. A sample of nearby pairs has been constructed from the UGC 
catalog (Nilson 1973) with the same projected physical separation; only a small percentage of the 
equivalent close pairs of galaxies are classified as peculiar (3 of 140) or uncertain (6 of 140). (This 
nearby pair sample is discussed further below.) 

4. IMPLICATIONS FOR MERGER RATES AND MASS EVOLUTION 

In this section we examine the consequences of our observations for simple merger models of 
galaxy evolution. We argue below that mergers must be taking place in our pair population, and 
from this infer that mergers must have taken place at a significant rate in the field population as 
well. Formally, none of our observations (mean redshift, colors, line strengths) demand that these 
mergers cause star formation. However, we argue below that our observations are in fact consistent 
with models in which mergers induce evolutionary brightening of both the paired galaxies and the 
field population. 

4.1 Pair Fraction and Merger Rate 

The fraction of galaxies that are members oi 6 < 6" pairs is about 14% of our sample (30 of 
196 for V < 21.5, and 52 of 410 for V < 22.5). We estimate that about 1/2 of these are physical 
pairs, based on our redshift data (3 of 7 physical pairs) and the measurement that cj{9) — 2 for 
these pairs. Therefore about 7% of the faint galaxies in our sample at z ~ 0.4 are in physical pairs 
with separations projecting to 19h~^ kpc. 

How does this result compare with the fraction of close pairs found by Burkey etal. (1994) in 
deep HST images? Their fraction of galaxies in pairs is 34%, for galaxies in the magnitude range 
18 < I < 22 and separations 0.5" - 4". If we assume uj{6) ~ 2 - 5 for 6* < 6" (as in Table 2), and 
extrapolate uj{9) to smaller separations using 6 = 0.8, we would predict fractions of galaxies in pairs 
around 12-14% - a factor of 2.6 x smaller than observed by Burkey etal. . Therefore the fraction 
of galaxies in close pairs, and the correlation function, is increasing much more steeply at small 
separations than a canonical power-law slope of 5 = 0.8 would predict. This is broadly consistent 
with the fact that our data also shows excess correlation power compared to a power-law fitted at 
larger separations. 

At separations of < 20h~^ kpc, interacting pairs should quickly merge. Toomre (1977) es- 
timates that the time to merge is typically 0.5 Gyr for such objects; this estimate should be 
relatively conservative, since the orbital time for a 200 kms""*^ rotation velocity is OAh~^ Gyr at 
this distance. The merger timescale for the population of faint galaxies as a whole is therefore 
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tmg{z = 0.4) ~ 0.5/~g Gyr, or 7.1 Gyr. For Hq = 50, this is essentially identical to the Hubble 
time at z = 0.4 (7.9 Gyr), implying that such a large rate of merging will have a significant impact 
on the average masses of galaxies. Note that all merger rates scale with the assumed 0.5 Gyr for 
pairwise mergers. This is likely to be a lower limit, with mergers occurring more slowly in groups 
where tidal fields help to delay merging. 

A low redshift sample with comparable pair properties to our faint galaxy sample can be 
constructed from the (angular diameter limited) UGC catalog, which is statistically complete to 
B = 14.5 (Nilson 1973). The average redshift (of the galaxies in the Nilson catalog which have 
redshifts) is 0.007 (2100 km s"^). At this redshift 19/i"^ kpc projects to 192". Restricting the 
catalog to objects with B < 14.5, with measured sizes in both B and R, we find 140 paired galaxies 
of 3066 in the sample, or a fraction of 4.6%. Taking all of these to be physical pairs, we estimate 
tmg{z = 0) = 11 Gyr, where we again assume that the timescale for an individual merger is 0.5 
Gyr (as for the faint pair sample). 

Although projection is not a problem for nearby galaxy pairs, some physically close galaxies 
will be moving at such high relative velocities that they are unlikely to merge. Of the 70 UGC 
pairs, there are 23 for which both galaxies have redshifts; only 13 of these have line of sight velocity 
differences less than 350 kms~^ and hence are likely to merge (the average velocity difference of 
the "< 350 kms"-*^ " group is 107 kms"-*^ , compared to 570 kms""*^ for the "> 350 kms"-*^ " group). 
Applying this fraction, 0.56, to the 140 close galaxies, reduces the number likely to merge to 79, 
which decreases f^g to 2.6%, and increases the merger timescale to 19 Gyr. It is interesting to note 
that 57% of the UGC close pairs are E or SO galaxies, whereas only 16% of the catalogue has these 
morphological types. Not surprisingly, 80% of the high velocity pairs are E or SO galaxies. It is 
noteworthy that none of the galaxies in the 13 close pairs with small pairwise velocities is classified 
peculiar, although there are 2 irregulars. 

The timescale for merging at the current epoch found here, 19 Gyr, is substantially shorter 
than implied in Toomre's (1977) discussion (but see Carlberg and Couchman 1989). The principal 
difference between these two results is that Toomre demanded that merging objects be severely 
distorted, with long tidal tails. In fact, there are many more interacting objects in the Arp and 
Madore (1986) catalog (by a factor of 5 or so) which are likely to merge, yet which lack such 
conspicuous evidence of interactions. 

4.2 Evolutionary Effects 

The similarities of the faint galaxy pair and field populations have two possible explanations. 
Either the presence of a close companion has no effect on the color or luminosity of a galaxy, or, if 
a close companion leads to interactions which boost the star formation rate (as argued by Larson 
and Tinsley 1978), then interactions must be so frequent that most faint galaxies in the field have 
recently undergone interactions. At redshift z ~ 0.3 the luminosity density, roughly proportional 
to the star formation rate, is about 3 times higher than at the current epoch (Broadhurst etal. 
1988, Lilly etal. 1991, Eales 1993). The following crude calculation shows that the small fraction 
of galaxies in close pairs, ~7% here, could be the source of all the enhanced star formation in the 
field. The Larson and Tinsley (1978) model for the B luminosity of a burst is approximately fit by 

Lk/M = 0.9t-°-^ (2) 
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where t is measured in units of Gyr, M/L is in solar units, and L is measured in the (rest-frame) 
B band. The ratio of the total luminosity of a galaxy having a star formation burst of age t}, to 
the average luminosity of post-starburst galaxies of maximum age t^ is 

(L) loi^ + Mo/M,' ^^ 

where T = tb/tm, and Mq/M^ is the ratio of the mass of the underlying pre-existing population to 
the mass of stars formed in the burst. 

Equation (3) allows an upper limit to be set on what fraction of an average pairs's stellar 
mass is being formed in the interaction, based on keeping Lb/{L) < 1.58 (O'T'S, the limit implied 
by the redshift distributions). For interactions occurring at a uniform rate, tb/tm can be set to 
the fraction of interacting galaxies, 7%. Using Equation (3) with tb/tm = 0.07 gives Mo/Mb > 12, 
which sets the upper limit to the fractional burst mass, Mb/MQ, as less than 8%. If the underlying 
old population continually forms stars, then it has a higher luminosity per unit mass than a old 
pure burst, which could easily lower this limit by a factor of 2. Therefore, induced star formation 
in close pairs cannot dramatically change the mass of a galaxy-close pairs at z ~ 0.4 are not 
"protogalaxies" . 

The cumulative effect of interactions over the lifetime of the galaxy could be a dominant 
source of its stellar mass, if each interaction does induce the allowed limit of star formation. That 
is, if 8% (or 4% for continuous star formation) of the stellar mass is formed in each interaction, 
and all the field galaxies are "post interaction" , which is allowed since the ratio of the interaction 
time of O.SGyr to the age of the universe at z ~ 0.4 is about 0.06, then the upper limit to the 
accumulated burst mass is 100 (50)%. Therefore, most of the mass of a typical galaxy could be 
built up as a result of interaction induced star formation (assuming that sufficient gas is always 
available). Greatly improved constraints will be possible with a larger sample. The data presented 
here includes the "interaction formation" hypothesis only as an upper limit. 

4.3 Redshift Dependence of Merger Rate 

The relative rates of merging at z = and z = 0.4 can be compared without knowing the 
timescale for merging. We estimate that 7. Oil. 5% of the galaxies at {z) ~ 0.4 are in pairs with 
the same physical separation and merger probability, and 2.6ib0.3% of the galaxies at z ~ 0. The 
ratio of merger rates at these two redshifts is therefore 2.7 it 0.7, where the uncertainties have been 
computed by assigning N^'^ errors to the number of objects at each redshift. (Note that the error 
in the fraction of objects in pairs at z=0.4 may be larger than derived from N^'"^ errors, because 
of uncertainties in the fraction of pairs that are physical.) Assuming a (1 -|- z)"* increase in merger 
rate with redshift therefore implies that m = 2.95 it 0.8, in good agreement with m = 2.7 derived 
by Burkey etal. (1994) from HST images. 

The merger rate index can be tentatively linked to Q; low fl gives very little growth of structure 
at late times, although m > 2 for any reasonable fi. The key assumption is that galaxies moving 
at less than some critical pairwise velocity will merge (Aarseth and Fall 1980); the critical pairwise 
velocity is comparable to the escape speed. Under the assumption that the pairwise velocity of 
encounter of two galaxies evolves as expected in hierarchical clustering, (1 -|- z)™", and that the 
critical velocity remains constant (not completely true for an evolving mass distribution) then it 
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is straightforward to show that the merger rate index, m, should vary as m ~ 4.5^^''^'^ (Carlberg 
1990) for an unevolving mass spectrum. The rate of change will be reduced if the galaxies become 
significantly lighter in the past. Allowing for a 32% mass difference between redshift and 0.4 (see 
below) changes the relation to m = 3.20.^'^'^ . For our measured m, the implied $7 = 0.87 it 0.4. The 
large random error means that this is not a very useful constraint on fi; clearly a larger sample of 
pairs could rectify this problem. The other merger model parameters have relatively little effect 
over this range of redshifts. Certainly the rapid rate of change in merger rate found here argues 
against a low Q Universe; in such Universes the growth of perturbations "freezes out" at high 
redshift (z ~ f^""*^), so that the present day merger rate is slower and depends more weakly on z. 



Figure 1: The merger rate - redshift relation. The z = 0.007 data comes from close 
pairs in the UGC catalog. The z = 0.4 results are from this survey. The lines use the 
theory discussed in the text, which is based on a generalization of the Press-Schechter 
formula. The size of the ellipses indicates the approximate spread in the data along 
each axis. The theoretical predictions (solid lines) depend primarily on tmg, which is 
assigned values of 13, 26, and 45 Gyr (top to bottom). The two closely-separated lines 
in the middle show the result of changing Vmg from 280 to 140 kms~^ . Integration 
under either of the middle lines implies that the average galaxy at z = 0.4 has a mass 
that is 68% of that of present-day galaxies. Extrapolating to z = 1 reduces galaxy 
masses to 36% of present-day values, assuming that low relative velocity pairs merge. 

The merger rates estimated from the UGC catalog and our survey are displayed in Figure 1, 
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along with several parameterizations of the theoretical merger rate-redshift relation (Carlberg 1990). 
The merger rate per galaxy is given as 

where 7(0, x) is the incomplete gamma function (the integral of a Maxwellian speed distribution), 

Vmg is the critical velocity for merging, taken to be 280 kms~^ , and Vp is the pairwise velocity 

dispersion of galaxies, taken to be 600(1 + z)*-"""'^-'' '^^"^^^kms""'^ (Carlberg 1990). The constant c„ 

normalizes the expression to t~g at t = to- We take n = —1.5 as is indicated by the correlation 

function and the CDM spectrum for mass scales of 3-10 galaxy masses. This value of n predicts 

that the pairwise velocity dispersion evolves as Vp{z) oc (1 + z)~^'^, a 25% drop out to redshift 0.4. 

The critical velocity for merging will evolve as the galaxy masses change; however for v^g oc M^'^ 

the inferred drop of ~ 30% in the average galaxy mass (see below) is sufficiently small that changes 

in Vmg/vp are caused primarily by changes in Vp. For small values of the tmg this expression exceeds 

the rate of halo mergers, 

_ 2 1 

predicted at larger z. The Rmg plotted in Figure 1 is the minimum of the two expressions, with 
Equation (4 )dominating in the < ^ < 0.4 range where our data lies. The scaling with redshift 
of Equation (4) should be about right, but the normalization is subject to a lack of almost any 
knowledge of the details of merging of galaxy size objects in a cosmological context. 
The change in the numbers of galaxies as a result of merging is 



dn 
~dt 



^ - RmgUg (6) 



Under the assumption that the baryonic mass content of galaxies may be rearranged but not 
changed, Eq. 4 implies that the characteristic mass of galaxies decreases as 

M, {z) = M, (0) exp [- / R^g dt] (7) 

Jo 

with redshift. We make a smooth transition between the two rates by using an inverse quartic 
interpolation. The model line shown in Figure 1 implies that at z = 0.4, the characteristic mass of 
galaxies M^{z = 0.4) = 0.69M*(0). Extrapolating to z = I predicts that M*(z = 1) = 0.36M*(0), 
which is indeed a dramatic change in the population. (However, it should be noted that the 
uncertainty in this result is very large - nearly a factor of 2 ~ and depends on the accuracy of the 
theoretical model.) These estimates of mass changes can be taken as upper limits, since the time to 
merge may well be longer than the 0.5 Gyr assumed, and since some close pairs may fail to merge, 
as a result of strong gravitational tides in clusters and groups. 

4.4 Companion Galaxies? 

The excess small scale correlations discussed in §3 (about a factor of 2 or 3 more galaxies at 
9 < 6" than expected, a result which is more prominent in V than I) have at least two possible 
explanations consistent with our observed redshift distribution: either they represent a true density 
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excess, or they are due to a selective brightening only of intrinsically low luminosity companions. At 
low redshift the correlation function appears to be consistent with a power law on all scales (Lake 
and Tremaine 1980, Davis and Peebles 1983). At somewhat higher redshifts quasar absorption line 
studies suggest that galaxies may have associated substructure (Wolfe etal. 1986, Steidel 1993). A 
luminosity excess for close pairs is predicted in models in which the changes in the field population 
luminosity density at z ~ 0.3 are due to galaxy interactions and merging (Rocca-Volmerange and 
Guiderdoni 1990, Cowie etal. 1991, Broadhurst etal. 1992, Carlberg and Chariot 1992, Carlberg 
1992). For an unperturbed luminosity function of the form L~°'dliaL, with a ~ 1 (approximately 
the normal Schechter function value) a brightening of a factor of 2, as allowed by our limits on the 
redshift distributions, increases the numbers above some limiting L by at least a factor of 2. 

Our data is consistent with either a density excess or a companion luminosity enhancement. 
A large, photometrically accurate imaging and spectroscopic survey will be necessary to accurately 
define the mean properties of these populations, which have large intrinsic dispersions. At the 
optimal magnitude for pair surveys, y ~ 21 — 22, there are only 100 pairs per square degree, 
requiring a great deal of sky to be surveyed. 



5. CONCLUSIONS 

We have 53 redshifts of galaxies drawn from a V < 22.5 photometric sample of ~ 1000 
galaxies. Close pairs are defined as objects closer than 6" on the sky, for which we have 14 
redshifts. The pairs are compared to a sample of nearby pairs assembled from the UGC catalog 
with the same projected physical separation. The principal result of this study is that the mean 
redshift of the pair population is very similar to the mean redshift of the field galaxies at the same 
magnitude limit. This limits the differential luminosity evolution between the field and pairs to 
less than ~ 0.5 mag difference. The small magnitude difference implies that if an average pair has 
interaction induced star formation, it only adds 8% more stellar mass on the average. Integrated 
over the lifetime of the galaxy this could account for the entire stellar mass, but as an upper limit, 
given the observational constraints we have derived. 

A second result is that the rise of the fraction of galaxies in close physical pairs implies that 
the merger rate-redshift relation is (1 + z)"^-^^^^-^ over the interval < z < 0.4. This result is 
somewhat shallower than that of Zepf and Koo (1989); however it should be more reliable because 
uncertainties concerning the redshift distribution of the pairs have been removed. When combined 
with theoretical modeling of the expected merger rate-redshift relation, our observations have two 
important implications. If all of these galaxies do merge, then the average mass at z ~ 0.4 is 
reduced by 32%, and by a factor of ~ 2.7 at redshift 1 (the latter result depending significantly 
on the accuracy of the model used for the extrapolation). Secondly, the slope of the merger-rate 
redshift relation indicates O ~ 0.87ib0.4. The combination of model dependencies and large random 
error makes this density estimate less than compelling. 

There is an excess correlation at projected separations of < 20h~^ kpc. The strong correlation 
at r < 20h~^ kpc could be either a local bound population, or the result of a luminosity enhancement 
of intrinsically faint galaxies in close pairs. 

There are several limitations inherent in this study. A larger photometric and redshift sample 
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would allow a conclusive check of the color differences between field and pairs, which are expected 
to be in the range of 10%, because the field is so much bluer than observed locally. In this study 
there is no conclusive direct indicator that mergers have taken place, or will take place, other than 
the (difficult to quantify) distortion morphology; nor is there any direct indication that pairs have 
excess star formation rates. The relative rate of merging is established to within an accuracy of a 
factor of two, but the absolute rate is still poorly known, and the impact of merging on the mass 
function hinges upon an integral over the absolute rate. 
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